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Euclidean strong coupling expansion of the partition function is applied to lattice Yang-Mills the- 
ory at finite temperature, i.e. for lattices with a compactified temporal direction. The expansions 
have a finite radius of convergence and thus are valid only for )3 < Pc, where j3t denotes the near- 
est singularity of the free energy on the real axis. The accessible temperature range is thus the 
confined regime up to the deconfinement transition. We have calculated the first few orders of 
these expansions of the free energy density as well as the screening masses for the gauge groups 
SU(2) and SU(3). The resulting free energy series can be summed up and corresponds to a glue- 
ball gas of the lowest mass glueballs up to the calculated order. Our result can be used to fix 
the lower integration constant for Monte Carlo calculations of the thermodynamic pressure via 
the integral method, and shows from first principles that in the confined phase this constant is 
indeed exponentially small. Similarly, our results also explain the weak temperature dependence 
of glueball screening masses below Tc, as observed in Monte Carlo simulations. Possibilities and 
difficulties in extracting jS^ from the series are discussed. 
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1. Introduction 

Contrary to weak coupling expansions, strong coupling expansions are known to be convergent 
series with a finite radius of convergence. In the early days of lattice gauge theory they were used to 
get analytical results for some physical quantities of interest, such as glueball masses or the energy 
density of lattice Yang-Mills theories. These calculations were done at zero temperature, i.e. at 
infinite volume and temporal extent A^, of the lattice. 

Here we calculate such series expansions for the free energy density and screening masses 
with an infinite spatial volume and a compactified temporal extension A', of the lattice. In this way 
finite temperature effects are generated, giving us the opportunity to study the physical, temperature 
dependent free energy density in the confined phase. The physical deconfinement phase transition 
then corresponds to a finite convergence radius of the series, which one may try to estimate from 
the behaviour of the coefficients. 



2. Free energy density 



2.1 Cluster expansion 

The partition function of the lattice Yang-Mills theory is given by a functional integration of 
the exponentiated Wilson action over the corresponding SU(N) group space, 



Z = J DU exp 



Y^:^{TrU + TvU^-2N) 



(2.1) 



2N 



An expansion in the lattice coupling j8 by group characters Xr{U) and a cluster expansion yields 
the free energy density ^ 



f = -j-lnZ: 



-61nco(i8)-i £ a{C)Y\^{Xir. 



(2.2) 



c=[x;') ' 

where Q. = V - Nt the lattice volume and cq is the expansion coefficient of the trivial represen- 
tation, which has been factored out. The combinatorial factor a{C) is introduced via a moment- 
cumulant-formalism, and equals 1 for clusters C which consist of only one so-called polymer X,. 



The quantity in equation ( [2.2| ) is customarily called a free energy, even at zero physical tempera- 
ture, because the path integral corresponds to a partition function if one formally identifies j8 with 
1/r. Here we are interested in a physical temperature T = \/{aNt), realized by compactifying 
the temporal extension of the lattice. The physical free energy is then obtained by subtracting the 
formal (A'; = °°) free energy, which is analogous to a subtraction of the divergent vacuum energy in 
the continuum. Thus the physical free energy density reads 



/(M,m)=M,m)-/(oo,„). 

The contributing polymers Xi have to be objects with a closed surface, since 



j dUXriU) = 5,0. 



(2.3) 



(2.4) 
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Figure 1: Graph contributing to the lowest order fi{Nt,u) of the expansion of the physical free energy 
density at finite temperature. 

This means the group integration projects out the trivial representation at each link. To calculate 
the group integrals one uses the integration formula 

I dUXr{UV)Xr(WU-') = Xr(yW). (2.5) 

For a more detailed introduction to strong coupling calculations we refer to [|l|]. 

2.2 Results 

The graph contributing to the lowest order of the free energy density is a tube of length Nt 
with a cross-section of one single plaquette. The contribution of these tubes together with inner 
plaquettes is 

SU(2): fi{Nr,u) = --^M^iV'fl^', (2.6) 
SU(3): MNr,u) = [Z;^' +c^'] , (2.7) 

where u, v and w are the expansion parameters of the lowest dimensional representations of the 
corresponding gauge groups, 

SU(2): u = ^ + ff{P^) v=g + ^(r), 

SU(3): i, = A + ^(^2) v = ^ + ^(r) >^ = ^ + ^(i8')> (2-8) 

and we have used the abbreviations 

a = 1 + 3v — 4u^ , 
b = 1 — 3m — 6v + 8w, 

c = 1 + 3m + 6v + 8h'- 18m^ (2.9) 

Higher order contributions consist of such tubes with local decorations of additional plaquettes 
either in the fundamental or in higher representations. For the interesting cases SU(2) and SU(3), 
these contributions up to the calculated orders are 



SU(2): /(A^„m)= - Am^V 



1 + 12A^,.4 - ^-^N,u' + Un^ + ^A^, ) 



,(2.10) 
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SU(3): m,u) 



Nt 
3 

— u 

N, 



l + \2NtU^ + 42NtU^ 



l + nN,u'^ + 30NtU^ 



115343, 



2048 

17191 
256 



1095327 



2048 

N,u^ + 63NtU^ 



N,u' 



(2.11) 



which are valid only for Nt >5. For smaller A^, there are modifications of these formulae coming 
from polymers with cross-sections larger than one plaquette. The complete results for Nt =2 and 
3 in SU(2) are 



N=2: f{2,u)=--u' 



110 



58472 
405 



64897681 
65610 



Nt = 3: f{3,u) 



.12 



.2 .^4 37966 . 843898 

1 - 6m^ + 50m^ —u'' H TTr:^U'^ 



(2.12) 
(2.13) 



135 405 
2.3 Free energy density as a glueball gas 

Recognizing the first orders of the corresponding glueball masses (see [g and [^) for SU(2) 

. 2 98 4 20984 6 151496 g ^^^^^ 
m(A++) = -41nM + 2M2-— m4___„6____j,8^ ^2.14) 



m(E^ 



and SU(3) 

m(A++) 
m(£++) 

m(r+-) 

one can write 
SU{2) 

SU{3) 



26 

-41n M + 2u^ — —u* + 



27 



13036 
405 ' 



28052 
243 ' 



(2.15) 



. ^ 2 ^' ^ 4 297 . 858827 ^ 47641149 7 ^ ^ 

-41n M - 3m + 9m^ - — M^ - lu^ - —U^ + -^r-r--d' H TT^T^u' (2. 16) 

2 2 10240 71680 ' 

. 2 27 , _ 4 153 . 1104587 ^ 29577789 . 

-4 In M - 3m + 9m^ + llu' H H m\ (2.17) 

2 2 10240 71680 ' 



9 . 98 4 33 = 36771 ^ 117897 7 

-41n M + 3m + -M -u + —u ^tttt" H , m , 

2 4 4 1280 448 



f{Nt,u) 
f{Nt,u) 



1 


"^-m(A++)iV, 






1 


"g-m(A++)iV, 







(2.18) 



(2.19) 



(2.20) 



corresponding to a gas of non-interacting glueballs in a hadron-resonance-gas model [ffl], where 



e r . 



(2.21) 



This is a rather remarkable result. It allows to see from first principles that the pressure p = —f 
is exponentially small in the confined phase, and it explains the success of the hadron-resonance- 
gas model in reproducing the confined phase equation of state. Since the partition function is not 
directly measurable in Monte-Carlo simulations, the pressure is usually obtained by the integral 
method [||], where the expectation values of derivatives are computed and then integrated numeri- 
cally, 

P 



P_ 

2^4 



= Nf dp' [So -St]., 
lo Jpo 
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Figure 2: Left: Plot of the pressure density p vs. j3 (L + M = 4 Fades for SU(2) and Nt = 3). The plot range 
corresponds to the confined phase up to the critical coupling. Right: Plot of p vs. T /Tc for SU(3) from 
Monte Carlo data 



with So = 6Po and Sj = 3{Pt + Ps), where Pq denotes the plaquette expectation value on symmetric 
lattices and Pt^s are those of space-time and space-space plaquettes with Nt < Ns- The lower inte- 
gration limit is usually set to zero by hand, arguing with an exponentially small pressure in the low 
temperature regime. Our results now justify this assumption from first principles. 

2.4 Phase transition 

Physical phase transitions limit the radius of convergence on the real j3-axis, signalled by a 
singularity in the full free energy. We model the full function from the series coefficients by Pade 
approximants [L,M] with 

IT ^ - ^+aiu + ... + aLU^ 
6o + oi" + --- + om" 

and search for the zeroes of the denominator. The resulting L+M = 2,3,4 Pade tables for Nt = 2,3 
with the nearest real singularities are shown in table (|l|). 

Zeroes uq of the Pade approximant which are very close to a singularity often indicate that the 
singularity is superfluous and disappears as the full funtion is approached. Hence, removing the 



SU(2): Nt=2 SU(2): A^, = 3 



[L,M] 


Uc 




\Uc — Uq\ 




[L,M] 


Uc 




\Uc — Uq\ 


[1,2] 


0.4033 


1.8227 


0.0899 




[1,2] 


0.3467 


1.5133 


0.0219 


[0,3] 


0.4675 


2.2201 






[0,3] 


0.5009 


2.4538 




[2,2] 


0.5201 


2.5981 






[2,2] 


0.4623 


2.1853 


0.2388 


[1,3] 


0.4684 


2.2262 






[1,3] 


0.4347 


2.0098 


0.1373 


[0,4] 


0.4684 


2.2261 






[0,4] 


0.4617 


2.1820 





Table 1: Zeroes of the denominator (uc) and the numerator (uq) of the [L,M] Pade approximants and the 
corresponding value of j3c. 
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singularities with a nearby zero, we obtain estimates for the critical couplings, which are not far 
from the Monte Carlo results jS,. = 1 .8800(30) for Nt=2 and jS^ = 2. 1768(30) for A^, = 3 [g]. 

The L+M = 4 Pade's for Nt = 3, SU(2), are shown in figure The spread in the curves gives 
an estimate of the systematic error of the approximants at that order. The exponential suppression 
in the confined phase as well as the onset of the pressure upon approaching Tc is reproduced by the 
strong coupling series. 



3. Screening masses 



3.1 Zero temperature 



Screening masses are defined by the exponential decay of the spatial correlation of suitable 
operators. We used plaquette operators in our calculations. Temporarily assigning separate gauge 
couplings to all plaquettes, the correlator can be defined as |0] 



C{z) = {TrUpM TrUp,{z)) = N^j^^\nZ{^ ,^,^) 



(3.1) 



At zero temperature the exponential decay is the same as for correlations in the time direction, and 
thus determined by the gluball masses, the lowest of which may be extracted as 

m = -lim- In C(z). (3.2) 

Z. 

The leading order graphs for the strong coupling series are shown in figure (^). This leads to the 
lowest order contribution: 

C{z) = Au'^'- = A&-"'^\ (3.3) 
Thus the leading order for the screening mass is given by 

= -41nM(j8). (3.4) 

3.2 Finite temperature 

The graph contributing to the lowest order of the difference between the screening masses at 
zero and finite temperature is shown in figure (|3|). To lowest order the mass difference is 

A;n,(^,) = m,(A,) -m,(oo) (3.5) 

(3.6) 

= -iNtu"'''-^ (3.7) 

Thus one can see that the finite temperature effect on the screening mass is very small below Tc, as 
is also observed in Monte Carlo simulations (for references, see [^). 
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Figure 3: Graphs contributing to the lowest order of the expansion of the screening mass at vanishing and 
finite temperature. The correlation-plaquettes are painted black. 

4. Conclusions 

We performed explorative studies of strong coupling expansions at finite temperature. Our se- 
ries for the free energy density is to the lowest orders consistent with a free glueball gas. This result 
justifies the neglect of the lower integration constant in numerical calculations of the equation of 
state by the integral method from first principles. Moreover, it gives an explanation for the success 
of the hadron-resonance-gas model in reproducing lattice data in the confined phase. By extrap- 
olating the power series via Pade approximants and looking for the zeroes of the denominator, it 
is possible to get estimates for the critical value /3c of the deconfining phase transition, although 
higher order terms seem necessary in order to obtain some accuracy here. Finally, glueball screen- 
ing masses show a weak temperature dependence in the confined phase, consistent with what is 
found in numerical simulations. 
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